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Abstract 

This paper continues an earlier study of giant planet migration, examining the effect of planet 
mass and disc viscosity on the migration rate. We find that the migration rate of a gap-opening 
planet varies systematically with the planet's mass, as predicted in our earlier work. However, 
the variation with disc viscosity appears to be much weaker than expected. 



1 Introduction 

This work is a continuation of lEdgaJ l|2007h . which showed that Type 11 migration (the migration 
of a gap-opening planet) was not as simply as commonly thought. Most previous workers assumed 
that, once a planet opened a gap in a disc, it would migrate along with the general viscous evolution 
of the disc itself. One particular prediction of this approach is that the migration rate of the planet 
should only vary with the dis c viscosity. The d i sc mass should not affect the migration rate, unless 
th e disc m ass were very low (|Sver and Clarkj . Il995l : llvanov et al.l . Il999l . examined this limit). 

lEdgaij ((2007) demonstrated that the migration rate of a gap-opening planet was directly pro- 
portional to the disc surface density, in contradiction of the usual assumption. The planet is not 
a test particle within the disc, but must exchange angular momentum with it. As a result, the 
mi gration rate will co ntinue to depend on the various masses involved. 

Ilda and Liiil (|2004l ) suggested an alternative prescription for determining the Type II migration 
rate. They balanced the angular momentum change of the planet with the maximum viscous 
couple in the disc. Following this idea, lEdgaij (|2007t ) predicted that the Type II migration rate 
should be 



where the planet's semi-major axis is a, the local disc surface density is E, this gas viscosity is f 
and the planet's mass is A/p. In this paper, we shall study the variation in migration rate of a 
gap-opening planet with the planet's mass and the disc viscosity. 



2 Numerical Details 

Our basic technique is identical to that of lEdgaij l|2007h . using the Fargo code o f lMassetll|2000al lbl). 
Fargo is a simple 2d polar mesh code d e dicat ed to disc planet interactions. It is based upon 
a standard, ZEUS-like l|Stone and Normanl Il992h hydrodynamic solver, but owes its name to the 
Fargo algorithm upon which the azimuthal advection is based. This algorithm avoids the restric- 
tive timestep typically imposed by the rapidly rotating inner regions of the disc, by permitting 
each annulus of cells to rotate at its local Keplerian velocity and stitching the results together 
again at the end of the timestep. The mesh centre lies at the central star, so indirect terms 
coming from the planets and the disc are included in the potential calculation. We make use of a 
non-reflecting inner boundary, to prevent reflected waves from interfering with the calculations. 
The pitch angle of the wake is evaluated in the WKB approximation. The inner ring of active cells 
is then copied to the ghost cells, with an azimuthal shift appropriate to the pitch angle. Material 
which flows off the inner boundary is not added to the star (nor does the planet itself accrete). 
At the outer boundary, mass was added, to compensate for the viscous evolution of the system. 
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We use units normalised such that G — A/, + Mp = 1, while the planet's initial orbital radius 
is set at a = 1. References to times in terms of 'orbits' should be understood to mean "orbital 
times at the planet's initial radius." The grid extends between r — 0.4 and r = 2.5. Scaled to 
Jupiter's orbit, this grid roughly covers the area between the asteroid belt and Saturn's orbit. We 
assume a constant aspect ratio disc, with h/r — 0.05, and initially constant surface density, Eq. 
We set Eo through 

^ 2 

9disc — — l^j 

which provides a quick estimate of the disc's mass within the planet's orbit. We set g^isc = 4x 10~^, 
which ensures that our disc is always significantly more massive than the planets we consider. 

The gravitational effect of the planet on the disc is smoothed at 0.6 of the disc thickness at 
the planet's orbital radius: 

= -^^ (3) 

where e — 0.6h. There are two motivations for this, the first being the desire to avoid having a 
singularity wandering around the grid. The second is physical. The 2d approximation becomes 
poor close to the planet, where the vertical distribution of material becomes important. The actual 
distance of material from the planet ceases to be the well approximated by the in-plane distance, 
which would lead to the gravitational effect being over-estimated. Accordingly, we soften the 
potential over distances comparable to the disc scale height. However, this softening length is still 
substantially smaller than the expected gap size and the planet's Hill sphere. When calculating 
the torque the disc exerts on the planet, material from within the Hill sphere is subject to an 
exponential cut off, for similar reasons. At the start of each run, the planet is introduced gradually 
(over about one orbit), and is not initially permitted to migrate. This is done to minimise the 
effect of transients caused by the sudden appearance of a planet in a smooth disc. We released 
the planet to migrate after fOOO orbits. The computational grid is covered by 128 radial and 384 
azimuthal cells (all uniformly spaced). 

So far as possible, t his setup mirrors that used in the comparison project presented by 
Ide Val-Borro et "all (|2006l V In that comparison, the Fargo code was seen to give similar results 
to other codes used to study the disc-planet interaction problem. 



3 Results 

We will now show the results of our runs. We considered planet masses of 1.0, 1.5, 2.0, 2.5 and 
3 Mj (recall that Mj/M© ^ 10"^) and constant gas viscosities of 2 x 10"®, 8 x 10~®, 9 x 10"'^ 
and 10~^ in our units. These values ensure that the planets are always gap-opening. 

In Figure [1] we show how the migration rates vary with the planet's mass in discs of differing 
viscosity. Although the planets start to migrate at t = on these plots, recall that they had 
already completed 1000 fixed orbits, giving time for the disc to adjust to the presence of the 
planet. We cut the y-axis at 0.6, since at that point the m = 2 ILR of the planet encounters the 
edge of the grid. The migration rate of the planet therefore becomes unreliable. 

The general pattern is clear: increasing the planet's mass decreases the migration rate. This 
happens for all values of the disc viscosity, in accordance with the prediction of Equation [T] 
However, the value of the scaling constant is not necessarily the same. Integrating Equation [T] 
we find that it predicts a migration rate approximately 25% slower than we see here. One should 
compare with Edgar (2007), where the initial surface density profile was varied. The variation 
altered the detailed migration rates, although by less than the variation in the disc mass (which 
was varied in that paper). 

We will now study how the migration rate varies with the disc viscosity. We show some 
sample results in Figure (2] The data are the same as those used in Figure [1] but plotted in 
different groupings. We only show the results for the 1 Mj and 3 Mj cases; the other planet 
masses gave similar results. 

We see immediately that the Equation [1] is not predicting the variation with disc viscosity 
shown by our numerical experiments. The variation is far weaker in our runs than predicted 
above. The difference cannot be due to the simplifications made in deriving Equation [1] since 
that should not affect the scaling with v. It is possible that the viscosities used (especially 
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Figure 1: Planetary migration in discs with constant viscosity and initially constant surface density. 
The lines are marked with the planet mass in each case. Each plot corresponds to a different value 
for the disc's viscosity: 10"^ (TL), 9 x 10~*^ (TR), 8 x 10"*^ (BL) and 2 x 10"^ (BR). 
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Figure 2: Planetary migration in discs with constant viscosity and initially constant surface density. 
The lines are marked with the disc's viscosity. Top plot is for a 1 Mj planet, the bottom for a 3 Mj 
planet. These are the same data plotted in Figure [T] 
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V = 2 X 10"") are close to the intrinsic numerical dissipatiorQ of Fargo. However, even the 
difference between the v = 10~^ and = 9 x 10~® cases seems to be less than that predicted. 
Unfortunately, attempting to run with viscosities higher than 10~^ was not reliably possible with 
Fargo. 

4 Conclusion 

In this paper, we have continued the work of lEdgaJ l|2007l ), examining the migration rate of gap- 
opening planets. We have shown that the migration rate is inversely proportional to the planet's 
mass, as previously predicted. This strongly suggests that Type II migration should be regarded 
as an angular momentum exchange between the planet and disc, rather than the planet merely 
acting as a 'relay station' between the inner and outer discs. 

However, we have not seen the migration rate following the expected variation with disc 
viscosity. Our numerical experiments showed a much weaker variation than that predicted. It 
is possible that this is a numerical artifact; the intrinsic numerical dissipation of Fargo may be 
comparable to that produced by the physical viscosity we introduced. The best way of testing 
this would be to repeat our experimental set up using a different hydrodynamics code. We have 
no particular reason to believe that Fargo is misbehaving, but the possibility cannot be ruled 
out without an independent test. 

Type II migration is not as simple a process as often thought. The migration rate depends 
strongly on the the masses of both the disc and the planet, and this should be accounted for in 
monte carlo models of planetary system formation. 
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^We do not call this 'numerical viscosity', since there is no reason to believe that the dissipation acts exactly like a 
physical viscosity 
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